Using Superconducting Qubit Circuits to Engineer Exotic Lattice Systems 
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We propose an architecture based on superconducting qubits and resonators for the implemen- 
tation of a variety of exotic lattice systems, such as spin and Hubbard models in higher or fractal 
dimensions and higher-genus topologies. Spin systems are realized naturally using qubits, while 
superconducting resonators can be used for the realization of Bose-Hubbard models. Fundamen- 
tal requirements for these designs, such as controllable interactions between arbitrary qubit pairs, 
have recently been implemented in the laboratory, rendering our proposals feasible with current 
technology. 
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I. INTRODUCTION 

The idea of quantum simulation has led to a large num- 
ber of theoretical proposals and some remarkable exper- 
imental results over the past decade pj. For example, 
the transition between a supernuid and a Mott insulator, 
which is a classic condensed-matter-physics paradigm, 
has been investigated_in a controlled fashion using a gas 
of ultracold atoms 0- Based on the success of such 
experiments, there have been intense efforts to devise 
methods for the simulation of various physical problems 
whose theoretical analysis is challenging. One area that 
has been studied substantially in the theoretical litera- 
ture, with a number of fundamental questions still unan- 
swered, is that of lattice systems in arbitrary dimensions 
and topologies. Implementing such systems using nat- 
urally occurring systems in two or three dimensions is 
challenging because the required connectivity is incom- 
patible with the geometry of the physical space in which 
the spins or lattice sites reside. Here we propose the 
implementation of such exotic systems using electronic 
nanocircuits based on superconducting qubits (SQs) and 
superconducting resonators [3] . 

The basic elements required for the implementation 
of our proposals have all been demonstrated experimen- 
tally. In particular, in a recent experiment [||, Harris 
et al. fabricated a SQ circuit with N > 100 qubits 
and demonstrated fully controllable interactions in blocks 
of 8 qubits, i.e. essentially all pairs of qubits were cou- 
pled to each other with individually controllable coupling 
strengths. This ability to design at will tunable couplings 
between any pair of qubits is a crucial ingredient in our 
proposal for experimentally engineering exotic quantum 
architectures, which are hard to study otherwise. We in- 
vestigate a range of systems that can be implemented us- 
ing such SQ circuits with fully controllable interactions. 
The experiments that we propose should be realizable in 
the near future, particularly those related to the ther- 
modynamic properties of classical spin systems, where 
many-qubit coherence is not required. 

We start by reviewing the present-day technology of 



superconducting qubits, resonators and tunable couplers. 
We then go through our list of proposals. In partic- 
ular, we examine complex quantum systems in higher 
and nonconventional dimensions, such as spin lattices 
in nonintcger (fractal) dimensions [5|-|7j], and we con- 
sider the embedding of lattice systems on exotic topolo- 
gies, such as the Klein bottle and Mobius strip |8TllO|. 
Both types of systems (i.e. those involving nonconven- 
tional dimensions or topologies) could advance our un- 
derstanding of phase transitions [lll - [l3j and lead to prac- 
tical applications. For example, it is known that di- 
mensionality plays an important role in phase transi- 
tions [11| , and the ability to engineer any desired dimen- 
sion would provide a valuable experimental knob in their 
study. Furthermore, spin lattices in nontrivial topologies 
and in higher dimensions have been shown to be better 
suited for implementing passive quantum-information- 
protecting schemes than conventional systems [l5[ . Near 
the end of the paper, we discuss the implications of the 
in-situ tunability of the parameters and the possibility 
of performing novel topology- quench experiments. 



II. PHYSICAL IMPLEMENTATION 

Although various types of SQs could be used to im- 
plement the architectures that we propose here, we shall 
focus on flux qubits, since these were used in the ex- 
periments of Ref. [3] ■ The flux qubit is a superconduct- 
ing loop interrupted by a number of Josephson junctions, 
typically one to four [3|, as illustrated in Fig.QJa). When 
a magnetic flux close to half a flux quantum is applied 
through the loop, two quantum states (one with clockwise 
and the other with counterclockwise circulating current) 
are almost degenerate and form the qubit basis. These 
states are usually denoted by |t) and and the qubit 
can be thought of as a spin- 1/2 particle. The single-qubit 
Hamiltonian in the language of Pauli matrices a x,y,z is 
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FIG. 1: (Color online) Schematic diagrams of (a) a simple 
three-junction flux qubit, (b) a flux qubit with a tunable 
gap, (c) two flux qubits interacting through a tunable cou- 
pler, and (d) two resonators interacting through a tunable 
coupler. In (a) the external magnetic flux $ oxt threading the 
superconducting loop controls the parameter e in the single- 
qubit Hamiltonian (Eq. 1). In (b) the magnetic flux $ e con- 
trols the parameter e in the Hamiltonian, while $ a controls 
the parameter A. In (c) two flux qubits are coupled induc- 
tively to a common coupler, resulting in an effective coupling 
between the two qubits. The effective coupling strength J 
can be tuned through the magnetic flux $c. In (d) two LC 
circuits, i.e. resonators, are effectively coupled to each other 
through a tunable coupler. 



where A is the minimum gap at the half-flux bias point 
and e is the deviation from this point. It is worth noting 
here that e is an easily tunable parameter and that it has 
recently become possible to also tune A using the exter- 
nally applied magnetic fields 16], as shown in Fig. [ljb). 

If two qubits are placed next to each other, the 
magnetic-dipole interaction gives rise to a two-qubit cou- 
pling term of the form 



H[ n t.q— q — J&1&2: 



(2) 



where J is the coupling strength and the subscripts indi- 
cate the different qubits. For two directly coupled qubits, 
the coupling strength J is hxed by geometry and material 
properties. One could avoid this limitation and obtain an 
effectively tunable coupling strength by employing a cou- 
pler, an additional circuit element that mediates coupling 
between the two qubits [17| . A schematic diagram of this 
technique is shown in Fig. [He). By tuning the bias pa- 
rameters (e.g. the magnetic flux through the coupler's 
loop), one can effectively tune the inter-qubit coupling 
strength J. An additional advantage of using couplers 
is the flexibility allowed in their design, which leads to 
the ability to couple qubits that are separated by large 



distances and to produce coupling terms in any desired 
pairing of the qubits [3]. With the above architecture, 
one obtains the many-qubit Hamiltonian 



H, = 



(3) 



with at least the parameters and J,.j being tunable in 
situ. Further details concerning the circuit can be found 
in Ref. Q. 

While qubits are suited for the implementation of 
spin-lattice Hamiltonians, Bose-Hubbard Hamiltonians 
require harmonic-oscillator-like circuit elements, i.e. res- 
onators. These can be implemented either as lumped- 
element LC circuits, as illustrated in Fig. Hfd), or as 
coplanar- waveguide resonators Q. In both cases, the res- 
onator behaves as a linear oscillator with the Hamiltonian 



H osc — Tilo ( a) a + 



(4) 



where ui is the oscillator frequency and a) and a are the 
oscillator's creation and annihilation operators, respec- 
tively. Several experiments have demonstrated coher- 
ent coupling between resonators and qubits [3]. With 
the technology that has been developed in that context, 
there should be no difficulty in coupling resonators to 
each other using couplers, thus leading to tunable cou- 
pling with any desired pairing of the resonators. The 
interaction Hamiltonian is then given by 



int.osc— osc 



= ,/ 



CLjCli 



(5) 



where the subscripts indicate the different resonators. In 
writing this form for the Hamiltonian, we have assumed 
that J <C huj/n (with n being the typical number of ex- 
citations in each resonator) such that the rotating-wave 
approximation is valid. The resonators can now be seen 
as sites in a Hubbard-like model, and the number of ex- 
citations in any given resonator represents the number of 
(bosonic) particles occupying that site. One thus obtains 
the non-interacting Bose-Hubbard Hamiltonian 
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^2 J itj {a\aj + a)ai), (6) 

i,3 



with the parameters Jij being tunable in situ. The res- 
onator frequencies generally exhibit small deviations 
from the values specified when designing the circuit, and 
recent experiments have demonstrated resonators with 
tunable values of oj [1S| . 



III. NONCONVENTIONAL DIMENSIONS 

On the theoretical side, spin lattice systems in higher 
dimensions (d > 3) have been studied extensively in 
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the past [ll|. Interest in these higher-dimensional sys- 
tems stems from the important role that dimensional- 
ity plays in the physics of phase transitions and criti- 
cal phenomena, as well as in determining the magnetic 
and thermodynamic properties of materials. Over the 
years, a number of different theoretical methods have 
been used in studying higher-dimensional systems, in- 
cluding renormalization-group and Monte-Carlo meth- 
ods. However, the validity of these methods generally 
cannot be established rigorously 19]. Therefore, exper- 
imental investigation of higher-dimensional systems is 
highly desirable. The experimental realization of such 
systems is hindered, however, by the difficulties associ- 
ated with coupling spatially separated elements, such as 
distinct two-level systems. 

With SQ networks it should be possible to overcome 
the difficulty of arbitrary connectivity [20|. The Ising 
model in d — 1,2,3,... dimensions can be implemented 
using the connections illustrated in Fig. [2] In this con- 
text one might worry whether the crossing of some of the 
lines in Fig. 2 would be a problem. However, this is not 
the case, since such lines can be fabricated in different 
layers, similarly to what was done in Ref. Q. Naturally 
there would be physical limits on the number of layers in 
a realistic system, and one might worry that for larger 
numbers of qubits a larger number of layers will be re- 
quired. However, the number of layers does not depend 
on the system size, but only on the engineered effective 
dimension: in principle no overlapping connections are 
needed for d — 2, only two layers of couplers are needed 
for d = 3, and so on. 

Similarly to higher-dimensional systems, spin systems 
with nonintcger dimensions have also received much the- 
oretical interest over the past few decades but have not 
been implemented experimentally 0-0] • One example of 
interesting theoretical results in this context is that the 
Ising model can exhibit spontaneous magnetization at fi- 
nite temperatures on spin lattices with d < 2 Q ; another 
is the prediction that certain spin models change critical 
behavior from second to first order at specific nonintcger 
dimensions 0]. 

One way to obtain noninteger dimensions, without los- 
ing local structure, is to use fractal geometries. Indeed, 
the proposed SQ architectures can be used to implement 
the Ising model on a well-studied fractal that can be used 
to probe noninteger dimensions between d — 1 and d = 2, 
namely, the Sierpinski carpet [5|. The relevant construc- 
tion is illustrated and explained in Fig. [31 Note that 
a two-dimensional lattice of qubits without any qubits 
missing can be used to generate any desired dimension 
between 1 and 2. This is achieved by decoupling any un- 
needed qubits from the rest of the lattice using the tun- 
able couplers, effectively removing these qubits from the 
system and creating holes in their place. Alternatively, 
the Sierpinski carpet can be generated by engineering 
the connections between the different qubits on the chip, 
even if the qubits are not arranged in a two-dimensional 
lattice (as explained in Fig. 2). 




FIG. 2: (Color online) Engineering the effective dimensional- 
ity of a lattice system by changing the connections between 
qubits: (a) Nine qubits can be connected so as to form a 
linear chain [purple (black) connections] or a 3 x 3 square lat- 
tice [purple (black) and orange (gray) connections], (b) Eight 
qubits can be connected into a three-dimensional cube. 




FIG. 3: (Color online) First few generations of a Sierpinski 
carpet with dimension log [8]/ log[3] ~ 1.9: one starts from 
a square of 8 qubits, a 3 x 3 lattice with the middle qubit 
removed; then each site is replaced by the original, first- 
generation square; and so on. Any dimension between 1 and 
2 can be obtained by adjusting the total size and missing frac- 
tion in the first-generation square: if the length of the square 
is n qubits and the length of the missing core is / qubits, 
the effective dimension is log[n 2 — i 2 ]/log[n]. Starting from 
a three-dimensional cube, one can obtain any dimension be- 
tween 2 and 3, and so on. 



IV. EXOTIC TOPOLOGIES 

A substantial amount of theoretical work has been de- 
voted to the critical behavior of spin lattices with non- 
trivial topology ]8] . Experimental realizations have been 
possible with thin-film materials 0, [l(| , but with limited 
controllability of subsystems. An open question in this 
area is whether there is any general principle connecting 
the properties of a many-body system with the curva- 
ture or the topology of the surface on which the system 
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resides. It would therefore be desirable to have small- 
scale lattices with local addressability, whose underlying 
topology is also controllable. 

With the designs proposed here, this experimental goal 
becomes more feasible. For example, it would be possi- 
ble to engineer a spin system whose topology is either a 
torus or a Klein bottle. These two lattice topologies can 
be engineered as follows: we start from, say, the 3x3 
lattice shown in Fig.[2ja). In addition to the connections 
shown, we impose periodic boundary conditions along 
the vertical direction, i.e. we connect qubits (1-7), (2-8) 
and (3-9). The type of connections made along the hor- 
izontal direction then differentiates between a torus and 
a Klein bottle. For the torus we connect qubits (1-3), 
(4-6) and (7-9); while for the Klein bottle we connect (1- 
9), (4-6) and (3-7). Adding the connections required to 
implement either the torus or the Klein-bottle topology 
(or even both sets of connections together) requires the 
addition of only one layer of couplers, in principle. 

The simplest model one could study on these struc- 
tures is the Ising model, where nearest neighbors 
interact via a af <r| coupling term. One could also apply 
an effective transverse field, i.e. a af term in the Hamilto- 
nian, of strength A, on every qubit and obtain signatures 
of a quantum phase transition [l2| as the parameter A is 
varied. One of the most telling signs of the system pass- 
ing through a 'critical' value A c , is that the entanglement 
properties of the ground state wave function change quite 
drastically Even small systems of, say, N = 4, 9, 16 
qubits, will display cusps of increasing sharpness in some 
entanglement measure, such as the concurrence between 
two neighboring qubits. Such entanglement properties 
can be measured relatively easily in the proposed archi- 
tecture. Since the coupling strengths are all tunable, the 
qubits can all be decoupled from one another. Then, by 
selectively turning on and off certain coupling strengths 
and performing the required quantum gates, the neces- 
sary multi-qubit observables can be measured, e.g., fol- 
lowing quantum state tomography or entanglement wit- 
ness measurement protocols. Alternatively, if one could 
simultaneously couple all the qubits to a common res- 
onator, one could follow the ideas proposed in Ref. [l4| 



to detect the quantum phase transition through the re- 
sponse of the resonator to an external probe. 



V. BOSE-HUBBARD MODELS 

The architectures described above for spin lattice sys- 
tems in nonconventional dimensions and topologies can 
also be implemented using resonators, such that the re- 
sult is a Bose-Hubbard system. In such systems, one can 
study the transport properties of the excitations, which 
play the role of bosonic particles. In particular, one can 
investigate the superfluid-Mott-insulator phase transi- 
tion and the Anderson-localization phase transition. 

In order to investigate the superfluid-Mott-insulator 
phase transition, one needs inter-particle interactions. 
Such interactions can be engineered by coupling each res- 
onator to a qubit, similarly to the proposals of Ref. [2l| . 
In order to investigate the Anderson-localization phase 
transition, one needs disorder in the single-site energies 
or nearest-neighbor coupling strengths (i.e. hopping ma- 
trix elements). Engineering such disorder is straightfor- 
ward with superconducting circuits, where each site en- 
ergy and each coupling strength is individually tunable. 
A recent theoretical work also proposed the possibility 
of engineered time-reversal symmetry breaking in a sys- 
tem of coupled superconducting resonators [22]. This 
idea can also be implemented in the systems of arbitrary 
dimension or topology proposed in our work. 

An important quantity in the study of many-body 
physics is the two-point correlation function of the form 
{a\aj). The tunability of the coupling strengths in SQ 
systems enables one to measure this quantity relatively 
straightforwardly. One starts by turning all the couplings 
off. By measuring the number of excitations in each res- 
onator, and repeating the experiment a large number of 
times, one obtains (aja^). If then one couples only res- 
onators i and j with coupling strength J as described by 
H2, the number of excitations as a function of time is 
given by 



J 



(at (<)*(<)) 



( a l(0K(0))cos 2 (^ 



(at(0)a i (0)} S in 2 (^ 



2<4(0)ai(0) 



(7) 



One can therefore use the resulting oscillations in order 
to extract the imaginary part of (aj (0)aj(0)}. Repeating 
the same procedure with a it/ 2 phase shift applied to one 
of the resonators, allows the extraction of the real part 
of (4(0)^(0)). 



VI. QUENCH DYNAMICS AND 
TOPOLOGY-QUENCH EXPERIMENTS 



The in-situ tunability of the parameters in SQ cir- 
cuits enables one to perform quench-related experiments, 
where the parameters are changed from some initial con- 
figuration to a different one. What starts out being 
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the ground state or thermal-equilibrium state then be- 
comes an excited state that tends to relax to the new 
ground state or thermal-equilibrium state. How this re- 
laxation takes place, and whether it is possible at all 
under the constraints imposed by conservation laws has 
been a subject of extensive studies in the past [l3|, H3- 
The quintessential example of such quench problems is 
the Kibble-Zurek mechanism 24], which describes defect 
formation in systems that are quenched from one thermo- 
dynamic phase into another and is proposed as the mech- 
anism for pattern formation, such as galaxy creation, in 
the early universe. Using the SQ architecture discussed 
here, one can quench the parameters across any of the 
phase transitions mentioned above and analyze the re- 
sulting dynamics of the system. 

We also propose the idea of implementing topology- 
quench experiments, in which a system is initially em- 
bedded in one topology and, subsequently, its internal 
interactions are changed in order to obtain a different 
topology. The idea can be explained using the 3x3 lat- 
tice of Fig. [H[b), which can be turned into a torus or 
a Klein bottle, as explained above. The two different 
topologies differ by two of the connections implement- 
ing the boundary conditions: (1-3) and (7-9) versus (1- 
9) and (3-7). A topology quench can be performed by 
switching off the torus-generating connections and then 
switching on the Klein bottle-generating connections. It 
is known that the partition functions of the Ising model 
in these two different topologies are not the same [8| , and 
we would therefore expect to see signatures of the differ- 
ent orders as we change from one topology to the other. 
These may be manifested by the magnetization and ther- 
mal entanglement properties [I3j . In this case, it is the 
actual topology of the underlying lattice structure that 
changes as one set of interactions is turned off and an- 
other is turned on. The reordering of the quantum state 
via the Kibble-Zurek mechanism will consequently have 
to follow this change of topology of the lattice. 



VII. TYPICAL EXPERIMENTAL 
PARAMETERS 

Superconducting flux qubits typically have minimum 
gaps [Le. the parameter A in Eq. (1)] in the range 3-7 
GHz [3j. Recently, new qubit designs have made this 
parameter tunable in situ, with a minimum value of es- 
sentially zero [16j. The parameter A can therefore be 
tuned to any value between and 7 GHz. By changing 
the externally applied flux in the large qubit loop, the 
parameter e can be tuned to any value between and 
values higher than 20 GHz. 

The coupling strength between a qubit and a cou- 
pler can be designed through the kinetic inductance aris- 
ing from shared superconducting segments or Joseph- 
son junctions. Coupling strengths close to 1 GHz have 
been achieved pBj . and with straightforward parameter 
changes it should be possible to reach values of about 



2 GHz. Couplers typically have minimum gaps A in 
the range 10-20 GHz. These parameters can result in 
a maximum inter-qubit coupling strength on the order of 
a few hundred MHz [I?} (Note that, by the design of the 
circuit under consideration, the effective inter-qubit cou- 
pling strength is tunable). It is therefore possible to ex- 
plore parameter combinations extending from the regime 
where the single-spin energies (including both a x and a z 
components) are larger than the interaction energy to 
the regime where the opposite is true. This is the region 
where phase transitions are expected to occur, suggesting 
that the experimental investigation of these phase transi- 
tions and related critical phenomena in superconducting 
qubit circuits is feasible. 

Superconducting resonators with frequencies in the 
few-Gigahertz range, which is the natural range to use 
in qubit circuits, can be designed and fabricated with 
high controllability. In qubit-resonator circuits, cou- 
pling strengths in the range of 10-100 MHz are common, 
and recent experiments have achieved coupling strengths 
close to 1 GHz [26| . Since superconducting couplers have 
similar structure to qubits, similar resonator-coupler cou- 
pling strengths can be expected. The above numbers in- 
dicate that inter-resonator coupling strengths (i.e. hop- 
ping coefficients) and on-site inter-particle interaction 
coefficients on the order of a few hundred MHz should 
be achievable. Since the competition between inter-site 
hopping and on-site interactions governs the Supcrfluid- 
Mott-Insulator transition, one should be able to access 
the Superfluid-Mott-Insulator transition in the proposed 
architecture. Disorder of magnitudes smaller than or 
larger than hundreds of MHz can also be achieved, im- 
plying that the Anderson-localization transition can also 
be investigated. 

Typical decay rates for both qubits and resonators are 
on the order of 1 microsecond in superconducting cir- 
cuits. It has been a challenge to experimentally fabricate 
multiple qubits on one chip where all of the qubits have 
coherence times at that scale. However, it is expected 
that in the future long coherence in multi-qubit circuits 
will be possible. When that goal is achieved, it will mean 
that the overall coherence timescale will be long com- 
pared to typical parameters in the Hamiltonian, which 
are on the order of tens of nanoseconds or shorter. 



VIII. CONCLUSIONS AND OUTLOOK 

Facilitated by high levels of controllability and steadily 
improving coherence properties, superconducting qubits 
and resonators are finding various potential applica- 
tions in quantum information processing and condensed- 
matter physics. For example, there have recently been 
a number of proposals for using them as quantum sim- 
ulators [l|. In this context, the single qubit /resonator 
controllability and readout is a key advantage of super- 
conducting qubits compared to microscopic simulators, 
such as natural atoms or ions. 
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In this work, we have added to the list of feasible poten- 
tial applications of superconducting circuits the engineer- 
ing of lattice systems in arbitrary dimension and topol- 
ogy. In particular, we have proposed to engineer spin lat- 
tices in integer dimensions d > 3, fractal dimensions and 
nonconventional topologies. We have also discussed how 
Bose-Hubbard lattice systems in similar exotic dimen- 
sions and topologies can be implemented using supercon- 
ducting resonators. An advantage of SQ systems in this 
context, which is particularly exploited in our propos- 
als, is the high level of connectivity [27j between super- 
conducting qubits, resonators and hybrid qubit-resonator 
systems. Furthermore, the in-situ tunability of the pa- 
rameters allows for the design of quench, or even topol- 
ogy quench, experiments: in such experiments the inter- 
nal reordering of the system could be observed as the 
connections between lattice sites are changed. In this 
case, the very topology of the underlying lattice could be 
quenched, thereby opening the way to a rather different 
type of quantum quench experiments. 



We should emphasize that our proposal for investigat- 
ing signatures of phase transitions in the Ising model does 
not require multi-qubit quantum coherence. Addition- 
ally, fluctuations in the parameters can be tolerated for 
purposes of analyzing the presence or absence of phase 
transitions. Such experiments should therefore be eas- 
ier to realize in the near future. Experiments combining 
both scalability and long multi-qubit coherence times are 
expected in the coming few years, at which point the in- 
vestigation of quantum phase transitions and critical phe- 
nomena using superconducting lattice systems can also 
be realized. 
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